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ABSTRACT 

Archival X-ray spectra of the four prominent single, non-magnetic O stars 
C Pup, C Ori, ^ Per, C Oph, obtained in high resolution with Chandra HETGS/MEG have 
been studied. The resolved X-ray emission line profiles provide information about the 
shocked, hot gas which emits the X-radiation, and about the bulk of comparably cool stellar 
wind material which partly absorbs this radiation. In the present paper we synthesize X-ray 
line profiles with a model of a clumpy stellar wind. We find that the geometrical shape of the 
wind inhomogeneities is important: better agreement with the observations can be achieved 
with radially compressed clumps than with spherical clumps. The parameters of the model, 
i.e. chemical abundances, stellar radius, mass-loss rate and terminal wind velocity, are taken 
from existing analyses of UV and optical spectra of the program stars. On this basis we also 
calculate the continuum absorption coefficient of the cool-wind material, using the Potsdam 
Wolf-Rayet (PoWR) model atmosphere code. The radial location of X-ray emitting gas 
is restricted from analyzing the fir line ratios of helium-like ions. The only remaining free 
parameter of our model is the typical distance between the clumps; here we assume that at 
any point in the wind there is one clump passing by per one dynamical timescale of the wind. 
The total emission in a model line is scaled to the observation. There is a good agreement 
between synthetic and observed line profiles. We conclude that the X-ray emission line 
profiles in O stars can be explained by hot plasma embedded in a cool wind which is highly 
clumped in the form of radially compressed shell fragments. 

Key words: X-rays:stars, starsiindividual: ( Pup, ( Ori, ^ Per, Oph 



1 INTRODUCTION 

Young and massive 0-type stars possess strong stellar winds. The 
winds are fast, with typical velocities up to 2500 km/s, and dense, 
with mass-loss rates M>10~'' Mq /yr. The driving mechanism for 
the mass-loss from OB stars has been identified with r adiation pres- 
sure on spectral lines jCastor. Abbot & Kleinlll975h . However, it 
was pointed out early fLucy & Solomonl ll97(]|) and later further 
investigated l^wocki & Rvbicki 1984) that the stationary solution 
for a line-driven wind is unstable; small perturbation grow quickly 
and result in strong shocks giving rise to the production of hot gas 
and emitting of X-rays. 

X-rays from hot OB type stars were discovered by the Ein- 
stein X-ray observatory ^Harnden et ak 1979; Seward et al. 197^. 
A large number of OB type stars was detected by the Rosat All 
Sky Survey (Berghofer, Schmitt & Cassinelli 1996), providing firm 
evidence that the X-ray emission is intrinsic for O star winds. The 
low resoluti on Rosat spec trum of the 0-type supergiant Pup was 
modeled bv iHillier et alj 1 19931) . In orde r to explain the app arent 
weakness of absorption in the spectrum, iHillier et alj il993h pos- 
tulated that the X-rays are emitted far out in the stellar wind, at 
distances larger that 100 J?*, by an optically thin shell expanding 
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at constant velocity. The emission from such shell would produce a 
broad rectangular box-like line. 

A box-like line profile is distinctly different from a skewed 
triangular line expected from radiation a rising in the wind accel - 
eration zone a nd attenuated in the wind iMacFarlane et aljll99lh . 
Ilenacd i200ll) demonstrated that the asymmetry of the line corre- 
lates with the wind opacity and density. In a more optically thick 
wind the lines are more blue-shifted. Since atomic opacity is wave- 
length dependent, the lines are expected to differ in shape across 
the spectrum, with larger blue-shift at wavelengths where the wind 
opacity is higher. 

X-ray generation in the wind acceleration zone is predicted 
in the model by Feldmeier, Puis & Pauldrach (1997), who showed 
from hydrodynamic simulations that dense cool shells of gas form 
in deep wind regions. Collisions of these shells can lead to strong 
shocks. X-rays originate from radiatively cooling zones behind the 
shock fronts. The hot X-ray emitting plasma is thermal, has a small 
filling factor, and is optically thin. The produced X-rays can be 
significantly attenuated by the overlying cool stellar wind. Hence 
the blue-shifted lines are expected. 

High-resolution X-ray spectra of O-type stars were ob- 
tained for the first time by the Chandra X-ray observatory 
iWaldron & CassinellilEoOll) . The resolved emission line profiles 
offered a stringent test for the theory. The analysis of the line ra- 
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tios from He-like ions in C, Ori and C, Pup constrained the regions 
of X-ray formation relatively close to the stellar surface. The line 
width was measured to be less than corresponding to the maximum 
velocity, confirming that the line emission originates in the wind 
acceleration zone . Line blue-shifts wer e also detected, indicating 
wind attenuation <Cassinelli et all200ll) . 

However, the data revealed that the line blue-shifts do not 
change significantly f or lines at different wavelengths in contrast 
to wh at was expected JWaldron & Cassinelll2o'oil : lcassinelli et alj 
I2OOII) . Moreover, the line fitting using the "standard model" of 
smooth homogeneous wind was not able to to reproduce the 
lines observed in non-magnetic single O stars unless with signif- 
icantly reduced m ass-loss rates (Kramer, Cohen & Owocki 2003, 
ICohen et all2006ft . 

There is a good reason to believe that mass-loss rates from 
O stars are need to be revised downwards. A plenitude of obser- 
vational and theoretical evide nce indicates the clum ped inhomoge- 
neous nature of O star winds iRepolust et alj2004l Bouret, Lanz & 
Hillier 2005, Fullerton, Massa & Prinja 2006). The empirical mass- 
loss rates of WR st ars are generally reduced by factors of a few 
in a clumped wind Mamann & Koesterkelll99% but it is still not 
known by how much exactly in case of O stars jpuls et alj200d) . 

Radiative transfer in a stochastic medi um is distinctly d iffer- 
ent from the transfer in a smooth medium iPomranin3ll99ll) . es- 
pecially when the clumps are not necessarily optically thin. There- 
fore, it is inconsistent to apply the smooth- wind formalism in order 
to model spectral lines in a wind where the mass-loss rate estimates 
are based on the assumption of wind clumping. 

The transport of X-rays in a clumped stellar wind was studied 
analytically and numerically by Feldmeier, Oskinova & Hamann 
(2003) and Oskinova, Feldmeier & Hamann (2004), respectively. 
Motivated by the prediction of hydrodynamic simulations that the 
cool-wind clumps are radially compressed, we considered the case 
of anisotropic wind opacit y, where the optical thic kness of a clump 
depends on its orientation. iFeldmeier et al.l i2003h showed and ex- 
plained that the assumption of flattened wind clumps leads to 
nearly symmetric and blueshifted li ne shapes. A correspon ding 2D 
stochastic wind model developed in IOskinovaet'ai]<2004 demon- 
strated that clumping effectively reduces the wind opacity and its 
wavelength dependence. The latter effe ct leads to a similarity of 
the shapes of all X-ray lines. Recentlv. lOwock i & Cohen (2006) 
used a similar model, but restricted to the specific case of isotropic 
clumps. In the latter model a symmetric line profile results only 
when the mean separation between optically thick clumps is quite 
large, of the order of a few stellar radii. 

The main objective of the present paper is to compare, for the 
first time, the model lines from clumped winds with the lines ob- 
served in the X-ray spectra of the prominent single non-magnetic 
O stars C, Pup, ^ Per, C, Ori, and C, Oph. 

Our basic strategy is to adopt the best available parameters for 
the stars and their winds, compute on their basis the model line 
profiles, and compare them with the observation. Thus we do not 
infer the model parameters from the line fitting, but instead answer 
the question whether the observed lines can be described in the 
framework of the X-ray transport in stochastic media. 

We utilize the most recent empirical mass-loss rates of our 
sample stars, and evaluate the mass-absorption coefficient in the 
wind using state-of-the-art non-LTE atmosphere code. Observed 
X-ray spectra are analyzed to infer the zone of X-ray emission 
and the wind velocity field. The only remaining free parameter of 
our model is the average time interval between subsequent clumps 



Table 1. Chandra ACIS-S HETGS observations of prominent O-type stars. 



Name 


ObsID 


log(AfH) 


Flux^ 


exposure 






[cm- 2] 


[erg/s/cm^ ] 


[ksec] 


CPup 


640 


20.00 


9.0 X 10-12 


67 


COri 


610 


20.48 


6.0 X 10-12 


59 


COri 


1524 


20.48 


6.0 X 10^12 


13 


^Per 


2450 


21.06 


1.2 X 10^12 


160 


COph 


4367 


20.78 


3.1 X 10-12 


48 


COph 


2571 


20.78 


3.1 X 10^12 


35 



1 Flux not coiTected for interstellar absorption 

passing the same point in the stellar wind. For this parameter we 
adopt the wind flow time, Rt/voo- 

With all model parameters defined, we model each line and 
compare it with the observation. 

The paper is organized as follows. Observational evidences 
of strong wind clumping are discussed in Section 2. The observa- 
tions and broad-band spectral properties of our sample stars are pre- 
sented in Section 3. The analysis of line ratios of He-like ions, con- 
ducted in order to estimate the location and temperature of the X- 
ray emitting plasma, is presented in Section 4. Section 5 describes 
the stellar atmosphere model, and Section 6 briefly summarizes our 
radiative transfer technique for X-rays in a clumped stellar wind. 
The comparison of modeled and observed line profiles is done in 
Section 7, and a discussion follows in Section 8. Conclusions are 
presented in Section 9. 



2 CLUMPED WINDS OF EARLY-TYPE STARS 

There is mounting observational evidence of strong inhomo- 
geneity of stellar winds. Stochastic variable structures in the 
He II 4686 A emission line in C, Pup were revealed by Eversberg, 
Lepine & Moffat (19 98), and explained a s an excess emission 
from the wind clumps. iMarkova et alJi20o3) investigated the line- 
profile variability of Ha for a large sample of O-type supergiants. 
They concluded that the properties of the Ha variability can be ex- 
plained by a wind model consisting of coherent or broken shells. 
iBouret et alj i2005h conducted a quantitative analysis of the far- 
ultraviolet spectrum of two Galactic O stars using the stellar atmo- 
sphere code CMFGEN. They have shown that homogeneous wind 
models could not match the observed profiles of O V and N IV and 
the phosphorus abundance. However, the clumped wind models 
match well all these lines and are consistent with the Ha data. This 
study provided strong evidence of wind clumping in O stars start- 
ing just above the sonic point. A reduction of mass-loss rates by a 
factor of at least 3 compared to the homogeneous wind model was 
suggested. A similar result was obtaine d from fitting stellar w ind 
profiles of the P V resonance doublet bv lFullerton et alj i2006h . A 
sample of 40 Galactic O stars was studied, and the conclusion was 
drawn that the mass-loss rates shall be reduced by up to one order 
of magnitude, depending on the actual fraction of P V in the wind. 
The discordance was attributed to the strong clumping of the stellar 
winds. 

X-ray observations of massive stars allow to probe the column 
density of the wind. An analysis of the X-ray emission from the 
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Figure 1. De-reddened HETGS/MEG±1 spectra of prominent O stars. 



O stars 3 Ori iMiller et alJl2002h and ( Pup jKramer et aDl2003h 
have shown that the attenuation by the stellar wind is significantly 
smaller than expected from standard homogeneous wind models. 

Similar conclusions are reached from the analysis of X-rays 
from colliding wind binaries (CWB). Such systems consist of two 
massive early-type stars. The copious X-ray emission is produced 
in the wind collision zone. At certain orbital phases the X-ray 
emission from the colliding wind zone travel towards the observer 
through the bulk of the stellar wind of one companion. Deriving the 
absorbing column density from X-ray spectroscopy constrains the 
mass-loss rates. An analysis of XMM- Newton observations of the 
massive binary 7^ Vel was presented bv lSchild et alj i2'004) . They 
showed that the observed attenuation of X-rays is much weaker 
than expected from s mooth stellar wind m odels. To reconcile the- 
ory with observation, ISchild et alj i2004h suggest that the wind is 
strongly clumped. Similar conclusions were reached from Chan- 
dra and RXTE observations of WR 140 jPollock et all2005l) . Alike 
7^^ Vel, the column density expected from the stellar atmosphere 
models that account for clumping in first approximation only is a 
factor of four higher than the column density inferred from the X- 
ray spectrum analysis. 

Spectacular evidence of wind clumping comes from X-ray 
spectroscopy of high-mass X-ray binaries (HMXB). In some of 
these systems a neutron star is on a close orbit deeply inside the 
stellar wind of an OB star. The X-ray emission with a power-law 



spectrum results from Bondi-Hoyle accretion of the stellar wind 
onto the NS. These X-rays photoionize the stellar wind. The result- 
ing X-ray spectrum shows a large variety of emission features, in- 
cluding lines from H-like and He-like ions and a number of fluores- 
cent emission lines. |[Sako et al. ( 200^ reviewed spectroscopic re- 
sults obtained by X-ray observatories for several wind-fed HMXBs. 
They conclude that the observed spectra can be explained only as 
originating in a clumped stellar wind w here cool dense clumps are 
embedded in rarefied photoionized gas. IVan der Meer et alj <200^ 
studied the X-ray light curve and spectra of 4U 1700-37. They 
showed that the feeding of the neutron star by a strongly clumped 
stellar wind is consistent with the observed stochastic variability. 

These observational findings appear to be consistent 
with the stellar wind theory. The hydrodynamic modeling of 
Feldmeier et al. (•19971) (see also a pseudo 2D simulation of 
Dessart & Owocki 200^ predicts that the stellar winds are strongly 
inhomogeneous starting from close to the core, and with large den- 
sity, velocity and temperature variations due to the de-shadowing 
instability. This instability leads to strong gas compression result- 
ing in dense cool shell fragments (clumps). The space between 
fragments is essentially void, but at the outer side of the dense 
shells exist extended gas reservoirs. Small gas cloudlets are ab- 
lated from these reservoirs, and accelerated to high speed by ra- 
diation pressure. Propagating through an almost perfect vacuum, 
they catch up with the next outer shell and ram into it. In this col- 



4 L. M. Oskinova, A. Feldmeier, W.-R. Hamann 

Table 2. Stellar wind parameters 



Name 


Sp.Type 










M 


clumping^ 


ref 






[Hq] 


[kK] 


[km/s] 


[hi-] 


[A/oyr-i] 






CPup 


041 


18.6 


39 


2250 


1.6 


4.2 X 10-6 


yes 


1 


COri 


09.5I 


31 


32 


2100 


2.9 


2.5 X 10-6 


no 


2 


5 Per 


07.5III 


24.2 


35 


2450 


1.9 


1.2 X 10-6 


yes 


1 


COph 


09V 


8.9 


32 


1550 


1.1 


<1.8 X lO-'^ 


no 


3 



"yes" - clumping was taking into account for mass-loss estimates; "no" — mass-loss is not corrected for the clumping 
n Puis et al. amd) . 2) Lamers & Leitherer (1993) . 3^ Repolust et al. (2004) 



lision, the gas parcels are heated and emit thermal X-rays. The X- 
ray emission ceases when the wind reaches its terminal speed. In 
contrast, the cool fragment s are maintained out to large distances. 
iRunacres & Owockil i2005f) studied the ID evolution of instability- 
generated structures in the winds and demonstrated that the winds 
are inhomogeneous out to distances of 1000 i?, . Thus theory pre- 
dicts the existence of two disjunctive structural wind components 
- hot gas parcels that emit X-rays, and compressed cool fragments 
that attenuate this radiation. 



3 HIGH-RESOLUTION X-RAY SPECTROSCOPY OF 
SINGLE O STARS 

A growing number of 0-type stars has been observed with the high- 
resolution spectrographs on board of Chandra and XMM-Newton. 
The purpose of this paper is to study the X-ray line profiles, there- 
fore we selected the brightest stars with well-resolved lines. The 
HETGS/HEG detector on board of Chandra has a spectral resolution 
AA = 0.012 A, while HETGS/MEG has AA = 0.024 A. However 
the effective area of HEG is smaller than that of MEG. The RGS on 
XMM-Newton has larger effective area but more modest spectral 
resolution. For consistency we concentrate here only on the first 
order co-added HETGS/MEG spectra. 

We consider only single stars because in a CWB the geom- 
etry of the system can affect the line shapes (Henley, Stevens & 
Pittard 2003). However, there is no way to make sure that an O 
star is definitely single. Runaway O stars are distinguished by an 
almost complete lack of multiplicity (Hoogerwerf, de Bruijne & 
Zeeuw, 2001). Therefore, we have selected all three single O-type 
runaways observed by Chandra HETGS - (" Pup, ^ Per, and C, Oph. 
These stars are relatively X-ray bright and observed with expo- 
sures long enough to accumulate high quality spectra as shown 
in Table Q and in Fig. Q We h ave also considered a fourth star, 
Ori, albeit it is a known binary jHummel et all200Ch . because the 
initial interpretation of X-ray emission line profiles from this star 
caus ed doubts in the validity of t he shock model of X-ray produc- 
tion <Waldron & Cassinelll200'll) . The stellar and wind parameters 
of o ur sample stars are compiled in Tabl e|2| 

IWaldron & Cassinelli <200lh and ICohen et alj <200d) pre- 
sented a detailed analysis of Chandra observations of ^ Ori. 
ICassin elli et aL ( ,2001) and Ki'amer et al. { 200^ analyzed Chandra 
spectra of C Pup. XMM-Ne wton observations of C Pup were an- 
alyzed bv lKahnetaI]j200lh and lOskinova et alj <2006l) . IWaldrotJ 
i2005h reported observations of C, Oph. Chandra observations of 
5 Per were not yet published to our knowledge. 



We retrieved the public archival data of these four stars and 
extracted the spectra using the latest version of the Chandra soft- 
ware and the calibration data base. The de-reddened MEG spectra 
of the stars are shown in Fig.[T](see also XM M-Newton spectra of 
O stars in Fig. 8 from lPaerels & KahrJ iiool ). The emission lines 
seen in all spectra are resolved, with broader lines seen in the stars 
with higher wind velocities. The spectra exhibit emission from H- 
like and He-like ions of low and intermediate Z elements. The Si 
and Mg lines are most prominent in the ^ Pup spectrum, while oxy- 
gen lines dominate the spectrum of C, Ori. The ratio of nitrogen to 
oxygen line flu xes is h igher in the runaway stars compared to the 
C, Ori spectrum. lKahn e t al. ( 2001) found that in C, Pup the emission 
measures derived from the nitrogen emission lines are at least one 
order of magnitude larger that those of carbon and oxygen. They 
conclude that the elemental abundance ratios of nitrogen to oxygen, 
as well as nitrogen to carbon are considerably higher than solar in- 
dicating that the material has undergone CNO processing. This is 
in accordance with scenarios where runaways have previously been 
members of close binary systems and undergone mass exchange. 

The spectral energy distribution (SED) of all stars can be fitted 
using the standard collisional plasma models, e.g. apec. The spec- 
tral fits indicate the presence of multi-temperatu re plasma wi t h tem- 
peratures in the range 0.2 keV<fcTy <0.7 keV. Ixahn et alJ iioOll) 
reported the detection of continuum in the XM M-Newton spec- 
trum o f ^ Pup, although their fits were inconclusive. lOskinova et all 
i200d) analyzed archival C, Pup data and concluded that the lines to 
continuum ratio is in agreement with the inferred plasma tempera- 
tures. 

The emission lines seen in the spectra of O stars cannot be 
fitted by means of the stationary plasma emission line models im- 
plemented in the fitting software. The standard models are not ade- 
quate for the fast moving stellar winds, and cannot predict expected 
line profiles. 



4 LINE RATIOS FOR HELIUM-LIKE IONS 

The He-like ions show characteristic ''fir triplets" of a forbidden 
(2), an intercombination {x + y) and a resonance (w) line. In OB 
type stars, the fir line ratios can be used to constrain the tempera- 
ture and location of the X-ray emitting plasma. As was shown by 
lOabriel & JordaiJ jl969h the ratios R{n^) and G(T^) are sensitive 
to the electron density and to the electron temperature: 

ii(ne,T,,d) = (1) 
x + y 
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Table 3. Ratios of fluxes^ in the components of He-like ions estimated by fitting observed spectral features. 



5 



Star 


NelX 


Mg 


XI 


Si XIII 




R 


G 


R 


G 


R 


G 


CPup 


0.35 ±0.12 


1.3 ±0.3 


0.36 ±0.09 


1.05 ±0.15 


0.9 ±0.2 


1.1 ±0.2 


COri 


0.3 ±0.2 


0.9 ±0.2 


0.9 ±0.4 


1.1 ±0.3 


1.9 ±1.6 


1.4 ±0.8 


5 Per 


0.1 ±0.1 


1.4 ±0.5 


0.19 ±0.15 


1.17 ±0.35 


2± 13 


1±3 


COph 


1.1 ± 1.1 


1.13 ±0.35 


0.3 ±0.1 


1.0 ±0.3 


0.9 ±0.8 


0.6 ±0.3 



^ Fluxes are from the triple-Gaussian fitting of MEG spectral lines 
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Figure 2. MgXI line observed in the spectrum of C^Ori, gPer, and COph 
(co-added MEG±1). Vertical dashed lines indicate rest-frame wavelength: 
. - intercombination line, - forbidden line. The 
figures are connected for the 
l200lh. The solid lines are 



Aw - resonant line. Ax 



rest frame wavelengths in thi s and all following 
radial velocities taken from iHoogerwerf et al J 
triple-Gaussians fits. 
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G=(z+(x+y))/w 

Figure 3. Observed and theoretical R and G ratios for the Mg XI fir triplet. 
Observed values are denoted by t he star nam e. The curves represent the 
theoretical B. and G from Porguet et alj l200ll) for the corresponding T,.ad 
and the radii of fine formation. Observed values of G are consistent with 
electron temperatures Tc < 3 MK. Observed values of R roughly constrain 
the line emission regions between ^1.5 and < 5 for all stars. 



z + {x + y) 



(2) 



A strong radiation field can lead to a significant depopulation of the 
upper level of the forbidden line via photo-excitation to the upper 
levels of the intercombination lines. This is analogous to the effect 
of electronic collisional excitation, thus mimicking a high density. 
For the characteristic densities of O star winds the photo-excitation 
is the dominant mechanism for depopulation. Since the radiation 
field dilutes with distance from the stellar surface, the ratio between 
forbidden and intercombination line provides information about the 
distance of the X-ray emitting plasma from the photosphere. How- 
ever, this is possible only when the radiative temperature Trad of 
the radiat ion field at the wave length of the depopulating transition 
is known. lPorguet et'al]i200lh performed improved calculations of 
R and G line ratios for plasmas in collisional equilibrium and tab- 
ulated them for a wide range of parameters, including T^ad- 

Lines of O VII, Ne IX, Mg XI, and Si XIII are observed in all 
our sample stars, however OVII is well resolved only in ( Oph. 
Therefore, we use in our analysis only He-like ions of Ne, Mg, and 
Si. In the sample stars with large wind velocities, ( Pup and ^ Per, 
the w,x + y, and z components strongly overlap. Hence flux es- 
timates in each separate component are ambiguous. The situation 
is best for Oph, the star with the lowest wind velocity, where all 
components are well separated. To obtain the flux in the w, x + y. 
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Figure 4. The EUV spectium of (, Pup, observed with FUSE (thin line), 
compared to a PoWR mode! spectrum (dotted). The resonance doublets of 
C I V, N V and O VI are well reproduced, as well as the "forest" of iron-group 
lines. The O VI doublet can only be fitted with models when assuming that 
a diffuse X-ray field causes additional ionisation. 

and z components for a given ion, we have fitted all three compo- 
nents simultaneously as Gaussians. Figure|2| shows an example of 
the observed and fitted lines of a He-like ion. The derived R and G 
ratios for all stars are listed in TableO 

The tables from lPorguet et allHoOlh were used to constrain 
the temperature and location of the X-ray emitting plasma. The ra- 
tio G is sensitive to Te, while it is almost insensitive to Tie, Trad 
and the dilution factor W . In our stars G is large. The derived elec- 
tron temperatures are lower than the temperatures corresponding to 
the ionisation potential for the given ion. The emission measures 
are highest for the ions with lowest ionisation potential. This is in 
accordance with the analysis of differential emission me asure in 
^ Pup and ^ Ori conducted bv lWoidowski & Schuld <2005h . 

For electron densities less than ~ 10^^ cm~^ and in the pres- 
ence of a strong UV field, the ratio ii (forbidden to intercombi- 
nation lines) is not sensitive to the density but depends instead 
on the electron temperature Te, radiative temperature Trad, and 
the dilution factor W . The electron temperature is constrained by 
G. The radiative temperature of the radiation field at the wave- 
length of interest (see Table 3 in |Porauetetal]j200lh ') can be deter- 
mined from stellar atmosphere models. We employ PoWR model 
to calculate the radiation temperature as function of wavelength 
for different depth points in the wind. The radiative temperature at 
Af 864 A (Si XIII) as < 90 per cent of the effective temperature. 
For the lighter ions Af 1033 A (Mg XI), Af 1270 A (Ne IX), the 
wavelength of interest is on the red side of the Lyman jump in the 
stellar SED, and the radiative temperature is higher, roughly equal 
to the effective temperature. 

With Trad and Tc constrained, the ratio R can be used to es- 
timate the dilution factor W and, therefore, the radial distance r 
from the stellar photosphere at which the X-rays are produced. The 
ratios R and G and corresponding error bars measured in Mg XI 
for our sample stars are shown in the R — G diagram in Fig.|3l 
Theoretical values of R and G, calculated bv lPorguet et alj J200l]) 
for Trad = 30 kK and Trad = 40 kK, and for different distances 
r, are also plotted in Fig.|3| From comparison between observed 
and theoretical values of R and G we conclude that the Mg XI line 
is emitted by gas located between 1.5 J?* and 5 R* in all stars. It 
should be noted, however, that the grids presented in lPorouet et alj 
i200lh are quite coarse. Moreover, the observed lines, especially in 
^ Per, are very noisy. 

We constructed similar diagrams also for NelX and Si XIII. 
The latter line is sufficiently well resolved only in C, Pup and C, Oph. 
Comparison of R — G diagrams for the different He-like ions in our 




Wavelength [A] 

Figure 5. The radius were the radial optical depth of the wind becomes 
unity in dependence on the wavelength in the Chandra HETGS/MEG range. 
The calculations were done using the PoWR stellar atmosphere code 
(see text) with stellar parameters from Tablel2l The prominent edge at 
A 21.5 A is due to oxygen. The vertical dashed lines correspond to the wave- 
lengths of the studied lines (as indicated). 

sample stars allows to conclude that the X-ray emission originates 
not closer to the stellar photosphere than ~ 1.5 i?* and not further 
than ~ 10 i?*, i.e. in the wind acceleration region for all stars of 
our sample. 



5 MASS-LOSS AND STELLAR ATMOSPHERE MODELS 
OF SAMPLE STARS 

Full non-LTE stellar atmosphere models of O star winds are nec- 
essary to infer p arameters of general, cool stellar wind. Recently, 
IPuIs et aP yOOg) presented an up-to-date review of the effects of 
clumping on M determinations in O star winds. They analyzed 
Hq, IR, and radio fluxes of C, Pup and ^ Per (within a sample of 
19 Galactic O stars) using non-LTE models which account for line- 
blanketing and clumping. The derived stellar parameters are listed 
in Tabled The mass-loss rates are identical for H^, IR, and ra- 
dio determinations and constitute upper limit s, since the clump ing 
properties of the wind are still unknown (see lPuls et alj l2006h for 
details). The stellar pa rameters for C, Oph listed in Table|2|are from 
iRepolust et alj J2004) . The mass-loss is based on the analysis of 
photospheric lines from H and He and shall be considered as an 
upper limit. This is also the case for parameters of Ori in Table|2| 
jLamersetaJl999h . 

The aim of this paper is to model the attenuation of X-rays in 
a clumped wind, where the photons are absorbed only in clumps 
and can escape freely between them. The total mass of the wind is 
conserved and distributed stochastically, in the form of clumps. The 
optical depth of clumps for the X-rays is determined by the wind 
density and the mass absorption coefficient. In our clumped model 
the total radial optical depth of the wind is absolutely the same as 
for a homogeneous wind of the same mass-loss rate. 

We consider continuum absorption of X-rays in clumped cool 
wind caused predominantly by bound-free and K-shell ionisation 
processes. Introducing Cartesian coordinates, the z axis points to- 
wards the observer and the impact parameter p stands perpendicular 
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to the latter. Alternatively, a point {p, z) is specified by spherical co- 
ordinates, i.e. the radius r — p'^ + and the angle ■& — cos^ ^ ^ 
between the line of sight and the radial vector. Hence ^ = z/r. It 
is assumed that the time-averaged stellar wind is spherically sym- 
metric and expands with a radius-dependent velocity 

v{r) ^ (^1 - , (3) 

where Voo is the terminal velocity, and ro is chosen such that 
v{r — i?*) — 0.01 Woo. R* is the stellar (photospheric) radius. 
For a given mass loss rate A'l, the continuity equation defines the 
density stratification via p(r) — M/(47rr^u(r)). 

The wind can be specified by its total radial optical depth, 

/•oo 

r* = / K,„p{r) dr , (4) 
J R, 

where is the mass absorption coefficient, which describes the 
continuum opacity of the bulk of "cool" wind material at the fre- 
quency of the considered X-ray line. 

In order to provide realistic values for this k^, we em- 

i iloy the Potsdam Wolf-Ra yet (PoWR) model atmosphere code 
Hamann & GrafenejEooj) . For each star we calculate a model 
with the parameters from Table |2| and the appropriate chemical 
composition from the references given in Table |2| For instance, 
enhanced nitrogen and depleted carbon and oxygen iLamers et alj 
1 1999ft are taken for the ^ Pup model. Figure|4|shows that our model 
spectrum of Pup agrees well with the observed EUV spectrum. 

The O VI resonance doublet (see Fig. |4j can only be repro- 
duced by models if there is a diffuse X-ray field causing some 
"super ionisation", as already recognized by Cassinelli & O lsoni 
il979l) . In our model we adopt thermal bremsstrahlung emission 
from a hot embedded plasma component of 3 x lO'' K and adjust 
the filling factor such that the level of emerging X-rays corresponds 
roughly to the observation. 

In fact, Ki, chiefly depends on the chemical composition, but 
not much on other details of the model, nor on the radial coordi- 
nate. Table|4|gives at the wavelength position of the studied X- 
ray emission lines, for each star of our sample. In order to visualize 
how opaque or transparent the stellar winds are for X-ray radia- 
tion, we plot at which radii the radial optical depth becomes unity 
(Fig-EJ, for the wavelength range under consideration. At a few 
stellar radii, i.e. in the acceleration zone of the stellar wind where 
X-ray producing shocks should be located, the wind of ( Pup is op- 
tically thick for all considered X-ray lines (except perhaps Si XIV). 
In ^ Per and Ori, only the softer of the X-ray lines should en- 
counter optically thick continuum absorption. The outer wind of 
( Oph is almost transparent for the X-ray line photons. 



6 TRANSFER OF X-RAYS IN A CLUMPED STELLAR 
WIND 

Wind inhomogeneity alters the radiative transfer significantly. We 
have studied the effects of wind fragmentation on X-ray line for- 
mation analytically in the limit of infinitely many stochastically 
distributed fragments iFeldmeier et aljEo03ft . and numerically for 
a finite number of fra gments within an accelerating stellar wind 
tOskinova et all Eoo4) . The latter work describes in detail the 
model and the code used in the present paper. Here we briefly sum- 
marize the description of fragmented-wind opacity in order to in- 
troduce the parameters that will be used in the next section to model 
observed X-ray lines. 
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Normalized frequency 

Figure 6. Solid lines: model line profiles of O VIII in the spectrum of ( Ori 
for different fragmentation frequency, ng , as indicated. Dashed hne: smooth 
wind model hne. The line emission region is assumed between l.SiJ, and 
9/?,, velocity law is with /3 = 1. Frequency is measured relative to the line 
center and in Doppler units referring to the terminal wind velocity Vac ■ Flux 
is normalized to the hne flux in the unabsorbed case. All model lines in this 
paper aie convolved with a Gaussian according to the HETGS/MEG spectral 
resolution. 
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-1.0 -0.5 0.0 0.5 1.0 

Normalized frequency 

Figure 7. The same as in Fig.|6]but lines are normalized to the same maxi- 
mum. The average time between subsequent clump passage, <T>= 1/no, 
is indicated in the upper right corner 

For a homogeneous wind, the optical depth Tii(p, z) between 
a point {p, z) and the observer is 

, , K,^M f°° dz' 

Jz r^v(r ) 

with r' — y'' p^ + z'"^. 

We assume that close to the stellar core at distances smaller 
than -Rmiii < 1.5-R* the wind is homogeneous and no X-ray pro- 
duction is taking place. The wind instabilities start to develop at the 
distances larger than i?min. These instabilities result in the heating 
of some small fraction of the material to X-ray emitting temper- 
atures. The general stellar wind remains cool and is compressed 
into shell fragments. While general wind fragmentation persists till 
large distances -Rmax ~ 300i?, , the X-ray emission is ceased at 
the distance -Rcmis > Sii* . We assume that the wind is homoge- 
neous at distances r> -Rmax and take into account the absorption of 
X-rays by this part of the wind as well. 
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The hot material (Tx ~ 1 MK) emitting X-rays is permeated 
with the cool (r„ ~ 10 kK) wind fragments (i.e. clumps), which 
attenuate the X-ray emission. Both the hot and the cool gas com- 
ponent move outwards with the same velocity v{r) and distributed 
stochastically. 

We assume for simplicity that all parcels of gas which emit 
X-rays have the same temperature Tx and, on average, the same 
mass. The line emission is powered by collisional excitation and 
therefore scales with the square of the density. We assume that the 
hot gas parcels expand according to the continuity equation when 
propagating in radial direction. Hence, the emission of each hot 
parcel is Lq / {r^v{r)), with a constant Lo determined by the line 
emissivity. 

The mass and the solid angle of each cool wind fragment are 
also preserved. In order that the time-averaged mass-flux M of 
shells resembles that of a stationary, homogeneous wind, the radial 
number density of shells must scale with l/v(r): 

n{r) = — - . (6) 
v{r) 

The parameter no = n(r)v{r) defines the number of frag- 
ments passing through some reference radius (e.g. entering or leav- 
ing the wind) per unit of time and is constant due to the mass 
conservation. It does not depend on the distance and is therefore 
convenient model parameter. We refer to no as the fragmentation 
frequency. The inverse quantity, <T>, is the average time interval 
between two consequent shells passing the same point in the stellar 
wind. When <r> is equal to the flow time, Th = _R*/«oo, the av- 
erage separation between clumps is 1 i?* when v{r) = Voc- In the 
wind acceleration region, where the X-rays are emitted, the radial 
separation between clumps is smaller. All model lines presented in 
this paper are computed assuming no — T^^ ^ . 

When no is specified, the time averaged number of fragments 
in radial direction is 

{N,) = no / (7) 

As was discussed in lOskinova et alJi2004) . the average radial 
optical depth of a fragment located at distance r is 



— rad 



47r 



1 
no ' 



(8) 



We assume here that the cool fragments are compressed ra- 
dially and therefore aligned (this model sometimes is referred to 
as "Venetian blinds"). The optical depth of a flattened fragment j 
along line-of-sight depends on its orientation 

f. = ^, (9) 



where 



(10) 



Thus, the optical depth of each fragment, and theref ore the effective 
wind o pacity, is angular dependent. As discussed in lOskinova et alJ 
i2004) . the optical depth along the line of sight between point (p, z) 
and the observer in a stellar wind which consists of aligned frag- 
ments is 

(1 - e-^^- ) Hr')\ ^ + Mp, w), (1 1) 



t{p,z) 



no 



v(r') 



where Zmax ~ \/ -Rmax — P^, and dz' = dr'//^. In the spe- 
cific case of isotropic wind opacity, which was considered in 



Table 4. Mass absorption coefficient [cm^ g ^] at the wavelengths of 
X-ray emission lines from stellar wind models, for the radial range between 
1.5i?, andSK, 



Line 


Wavelength [A] 


CPup 


COri 


g Per 


COph 


Si XIV 


6.18 


18 


34 


11 


11 


MgXII 


8.42 


22 


44 


23 


23 


NeX 


12.14 


59 


73 


55 


57 


FeXVII 


15.01 


94 


108 


95 


95 


O VIII 


18.97 


171 


177 


170 


172 


N VII 


24.78 


93 


80 


93 


94 



lOwocki & Coherj i2006l) . factor /i(r) in omitted in Eq. ^9) and 
Eq.GB- 

Let us compare the optical depth in a fragmented wind 
(Eg. ll H with the optical depth in a homogeneous wind (Eq.|5}. In 
the limiting case of optically thin fragments {fY^ ^1) the expo- 
nent under the integral in Eq. <I 1> can be expanded. Substituting the 
average optical depth of a fragment defined by Eq. j9j the optical 
depth in the thin-fragment limit is the same as in a homogeneous 
wind. 

The dependence of the model line on no is demonstrated in 
Figs.|^ One can see from Eq. jsj that when no is small the frag- 
ments are opaque for X-rays. In this case the radiation can escape 
only between fragments and the shape of the line profile does not 
change as long as the fragments remain opaque. With increasing 
no the number of optically thin fragments increases, changing the 
line profile. For sufficiently large ?io, when all clumps are optically 
thin, the line profile looks like from the smooth wind model. 

In the limit of optically thick fragments {fY'^ S> 1) the optical 
depth becomes 



r(p, z) = no 



dr' 
D(r') ' 



(12) 



where the contribution from rh(p, z(-Rmax)) is omitted for clarity. 

The dependence of optical depth on the mass absorption co- 
efficient and M has disappeared in Eq. <I2> . Thus, in a limiting 
case of a wind which consists only of opaque fragments, the opac- 
ity is grey and is determined by the fragmentation frequency no 
and the velocity. 

In our stochastic numeric wind model the fragments may have 
different optical depths, and are not necessarily opaque. As can be 
seen from Eq. <1 It even when the optical thickness of the fragments 
is ~ 1 the dependence of optical depth r(p, z) on is weaker, 
compared to the homogeneous case (Eq.|5j. 

Hence, albeit the opacities ac^ differ significantly depending 
on the wavelength (see Table|4}, the optical depths t(p, z) for dif- 
ferent lines remains similar (more so in the case of optically thick 
clumps). Therefore the profiles of X-ray emission lines are similar 
across the spectrum. 
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Figure 8. Strongest emission lines observed in CP^p (co-added MEG±1) and the model lines. The solid lines are the model with no = 1.7 X 
10~^ (<T>= l.Ghr), /3 = 0.9, and the radiation originating between 1.5-R* and 5-R*. The dashed line is for the homogeneous model of the same M. 



7 COMPARISON WITH OBSERVED X-RAY LINE 
PROFILES 

7.1 Model specification 

The model is specified by mass-loss rate M, velocity field param- 
eters i}oo and (3, mass-absorption coefficient k^, and fragmentation 
frequency no, as well as by the radii of the line emission zone. 

We use the most recent published mass-loss rates for our sam- 
ple stars. In case of Pup and ^ Per M is corrected for clumping, 
while for ( Ori and C, Oph the clumping is not accounted for. Ve- 
locity field parameters are from spectral analyses referenced in Ta- 
bleU The mass-absorption coefficients are determined using the 
sc PoWR code. For all our models we adopt the fragmentation fre- 
quency which is equal to the inverse wind flow time, no — Voo/R* ■ 
The radii of the line emission zone are constrained from our analy- 
sis of fir line ratios in Sect. 4. Thus, no free parameter remains and 
a model is completely specified. 



depth unity at 2R,. This is in agreement with the estimates based 
on the^r analysis. 

The clumped wind model predicts that the shape of the emis- 
sion lines is similar a cross the spec t rum. T his prediction is con- 
firmed observationallv. lKramer et ai]^2003^ note from the analysis 
of observed lines that the amount of absorption inferred for dif- 
ferent lines in ( Pup and the line shapes are not sensitive to wave- 
length. We emphasize that the model lines shown in Fig.|8|are not 
selected best-fit models. These lines are computed using indepen- 
dently derived parameters. Comparing model and observation, we 
conclude that the stochastic wind model is capable to reproduce the 
observed lines sufficiently good. 

Some of the lines shown in Fig.|8| are contaminated 
by blends. The Nvil l ine is blended with N VI (A 24.898) 
jPollock & RaassenI [2006h . Similarly, FeXVII line is possibly 
blended with Fe XIX (A 15.08). Emission from Ovil (A 15.17) 
is also likely to co ntribute in the far red wing of FeXVII 
JOskinova et all200d) . 



7.2 X-ray emission lines in C Pup, ^ Per, ( Ori, and C Oph 

^ Pup. Figurejs] shows the strongest emission lines observed in 
^Pup. As can be seen in Fig.|5| the wind of CPup becomes op- 
tically thin for X-rays at the Si XIV A6.18 line just above 2R„. 
When the wind is optically thin above the line emission region, 
the emergent line should be nearly symmetric relative to the rest 
wavelength Ao. However, the peak of the Si XIV line is shifted to 
the blue, quite similar to other lines in Fig.|8| This indicates that 
the line photons must have originated below the surface of optical 



^Per. The strongest X-ray emission lines observed in ^Per are 
shown in Fig.|9| The data are quite noisy, however, the lines look 
similar to each other regarding broadening and blue-shifts. For 
comparison, we have over-plotted on the FeXVII and OviI line 
both smooth and fragmented model line profiles. The difference 
between them is relatively small, because the clumps in ^ Per are 
not optically thick. It appears that the red wing of the FeXVII 
line is above the model. This could be due to blending with lines 
from higher iron ions, or it may indic ate that the fragm entation 
frequency is somewhat lower in ^ Per. IPuls et all i200d) suggest 
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Figure 9. Lines of Ne X, Fe XVII, and O VIII observed in ^ Per (co-added MEG±1) and the model lines. The solid lines are the model with no = 1.4 X 10 * s ^ 
((T) = 1.9hr), /3 = 0.9, and the radiation originating between 1.5-R, and 7iJ,. Dashed hnes are homogeneous wind model. 
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Figure 10. Lines of Ne X, Fe XVII, and O VIII observed in C, Ori (co-added MEGitl) and the model lines. The solid lines are the model with no = 2.7 X 
10^^ s^^ (<T>= 2.9 hr ), /3 = 0.8, and the radiation originating between l.lRt and 9-R». The dashed line is for the homogeneous model of the same M. 



a much higher clumping factor in the outer wind of ^ Per than in nearly identical with the line computed using smooth wind model 
CPup. withM « 1.3 X lO"''M0yr"^ 



Ori. Figure fTol shows the best-resolved lines in the Ori spec- 
trum. The lines have similar shapes, they are blue-shifted and 
broadened. The model lines shown in Fig. llOI are calculated us- 
ing mass-loss which is most certainly is overestimated by factor 
of ~few. Therefore, the model lines are slightly shifted to the blue 
compare to the observed lines. As shown in Fig.0 the reduction 
of fragmentation frequency by factor of 3 can shift the model line 
to the red by about ~ 0.1«oo. The reduction of empirical mass- 
loss rate of C,On will also lead t o the improvement of the model. 
Lines of C, Ori were analyzed bv lCohen et alj i2006h by means of 
smooth wind model. They have found from the best-fit model of 
O VIII that K,jKI / 'i-nvoo ~ 0.26. Substituting the mass-absorption 
coefficient from Tableland stellar parameters from Table|2| the 
corresponding mass-loss rate is M ~ 1.3 x 10~^ M(7) yr~^, whic h 
is 20 times lower than estimated in iLamers & Leithereil jl993h . 
On the other hand, the clumped wind model provides the same 
quality of fit assuming M only twice as low as in Table|2| and 
fragmentation frequency equal inverse flow time. This is illus- 
trated in Fig. ll II which shows the observed and the four model 
lines of O VIII in C^On. The line from clumped wind model, with 
no = Tf[^ = 9.1 X 10"^ and M = 1.2 x 10"'' Mq yr"^ is 



( Oph. The wind of ( Oph is almost transparent for X-rays. As 
can be seen in Fig.|5| at the wavelengths of the Mg XII, Ne X and 
FeXVII the wind is optically thin for the X-rays already above 
~ 1.5 i?,. These lines are shown in Fig. ll2l Due to the nearly 
complete lack of absorption in the wind the lines provide detailed 
information about the location and conditions of the hot plasma it- 
self. The lines are symmetric and quite narrow, with broadening 
up to the ±0.8 Voo- The slight asymmetry is due to occultation by 
the stellar core. The line shapes indicate that emission originates in 
the wind acceleration zone. The line shape is best reproduced (by 
eye estimate) when /3 = 1.5. As can be seen from the two models 
over-plotted on the observed FeXVII line, the difference between 
smooth and fragmented models is negligible for the thin wind of 

COph. 

At the wavelengths of the O VIII and the N VII line the wind 
of COph is thin only above ~ 2.5 i?,. These lines are shown in 
Fig-El Both these lines have a low signal-to-noise ratio, and it 
is difficult to assess the influence of wind absorption on line pro- 
files. The O VIII line appears asymmetric, with the blue wing be- 
ing weaker than the red one, and our spherically-symmetric model 
cannot reproduce it well. However, it is not clear whether this line 
asymmetry is significant. The N VII line is slightly shifted to the 
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Figure 12. Lines of MgXII, NeX, and FeXVII observed in COpli (co-added MEGitl) and the model lines. The solid lines are the model with no = 
2.5 X lO^** s^^ (<T>= l.lhi'),/3 = 1 . 5, and the radiation originating between 1.5 R, and 9 iJ* . The wind optical depth is smaller than unity for all these 
lines. Consequently, all lines can be equally well be reproduced by a smooth-wind model (shown by a dashed line). 
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Figure 11. Observed (histogram) and model hnes of O VIII in f Ori. Arrows 
indicate models used to compute the lines. For smooth wind models, only 
mass-loss rate, M [Mq yr~^] is shown, for clumped wind models, both 
mass-loss rate and the inverse fragmentation frequency, <T>= is 
shown. Line with M = 1.3 X 10~^ M p, yr^^ is the best-fit smooth wind 
model, as found in lCohen et alj|20o3) . 



blue. To demonstrate how sensitive is the emission line profile to 
the location of the X-ray emitting plasma, we show two model lines 
for O VIII and N VII in Fig. 1131 For each ion, the model lines differ 
only by the assumed minimum radii of X-ray emission. It appears 
that the OviII and the Nvil lines originate from slightly further 
out in the wind than the lines shown in Fig. ll2l at ~ 1.9 R,. 

The ( Oph wind opacity is small and X-rays produced in the 
wind can escape without suffering significant attenuation. This al- 
lows to assess the intrinsic production of X-rays. We have calcu- 
lated the ratios of the X-ray luminosity to the bolometric and me- 
chanical (Mv'^/2) luminosity of stellar wind in C^Oph and com- 
pared it with these ratios for other stars. The results are given in 
Tabled 

The ratio of X-ray to bolometric luminosity in all stars in Ta- 
ble |5) is slightly smaller than the canonical 10^^ known for OB 
stars (Berghofer et al. 199^. Perh aps, this is because the l arge sam- 
ple of bright OB stars analyzed bv lBerghofer et alJh99d) is biased 
towards binary stars, which tend to be more X-ray bright. 

As can be seen from Table |5| the ratio of the X-ray and the 
mechanical luminosity is highest in ( Oph compared to the other 
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Figure 13. Lines of O VIII and Nvil observed in COph (histograms, co- 
added MEGibl) and the model lines. The solid lines are from a model where 
radiation originates between 1 .9 R* and 9 , while dashed model lines are 
for the radiation from 1.5 i?» — 9 Rt . 



stars with larger wind attenuation. We estimate that about 0.01% of 
the kinetic energy of the stellar wind of Oph is spent on the X-ray 
generation. 



8 DISCUSSION 

The result of our analysis is that the stochastic wind model lines are 
in good agreement with the lines observed in stellar X-ray spectra. 
This confirms the clumped nature of O star winds. 

So far, modern non-LTE steUar atmosphere models include 
clumping only in the first approximation of optically thin clumps. 
For clumps which are not optically thin, a non-LTE treatment re- 
quires to follow the transport of radiation inside the clumps, which 
is presently an impossible task at least in full detail. 

But fortunately, the problem of X-ray formation in the winds 
of massive stars is much simpler, since we assume that the X-ray 
emission of optically thin plasma is decoupled from the continuum 
absorption in cool wind. Therefore, the full non-LTE treatment is 
not needed, and the solution of pure-absorption radiative transfer is 
sufficient. 

Optically thin clumping in the winds reduces the empirically 
determined mass-loss rates. It is not a priory obvious how the pres- 
ence of optically th ick clumps may affect the mass-loss estimates 
iBrown et aLll2004 . There is no observational evidence that the 
clumps in the wind are all optically thin. On the contrary, pho- 
tometric and spectral variability indicate that opaque clumps are 
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Table 5. Ratios of X-ray, bolometric and mechanical luminosities 



Name 


D 




Lhoi 


-^mcch 


ix/ibol 


-^x/ -^mcch 




[pc] 


[erg/s] 


[erg/s] 


[erg/s] 






CPup 


460 


2.3 X 1032 


2.8 X 1039 


6.7 X ICS'S 


8.2 X IQ-* 


3.4 X 10-5 


COri 


500 


2.1 X 10^2 


3.0 X 10^^ 


3.5 X 10^^ 


6.8 X 10"* 


5.8 X 10-5 


^ Per 


850 


1.4 X 10^2 
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Figure 14. Observed (histogram, co-added MEG±1) and modeled Ne X line 
in f Pup. From left to right model line profiles are calculated assuming 
M = 43, 8.6, 4.3, 2.1, 0.43 x lO-'^MQyr-^. CoiTesponding normal- 
ization coefficient for the model line maximum to the observed line maxi- 
mum is r; = 2.6 X 10", 16.0, 2.6, 0.85, 0.23 



present in stellar winds. The outstanding question of "true" stellar 
mass-loss rates in clumped win ds is being extensively investigated 
at present (e.g. lPuls et alj2006l) . 



8.1 On the dependence of line profiles on mass-loss rate 

In our modeling we adopted mass-loss rates which, when available, 
account for the wind clumping. Since the mass is conserved in the 
wind, the optical depth in each clumps scales directly with M, and 
inversely with fragmentation frequency (Eq.|8}. The effect of wind 
clumping on the lines profiles is most pro nounced when clump s 
are optically thick for the incident radiation jFeldmeier et alj2003h . 
Therefore, the reduction of empirical mass-loss rates means that in 
order for clumps to be optically thick the fragmentation frequency 
should be smaller 

In Fig. ll4l we model the Ne X line in ^ Pup for the same 
fragmentation frequency, but different mass-loss rates. The low- 
est value, M » .43 x 10-"" Mq yr-^, is an estimate by 
iFullerton et al.l i200d) . based on the assumption that all phospho- 
rus is in the Pv ground level (g(P^^) — 1). Adopting this value, 
the effect of wind absorption on the line profile is smaller than ob- 
servered, as can be seen from the smaller blueshift of the model 
line compared to the observed line in Fig. ll4l With increasing M 
the line maximum becomes more and more blueshifted. 

The effect of increasing absorption on the line flux, which can- 



not be seen from Fig. 15 because the profiles are scaled to the same 
maximum, is dramatic. Actually, the profile for the highest AI is 
IC times weaker than for the lowest M, provided that the same 
flux has been originally released in the X-ray emitting zone. Even 
for the optically thin Fe XVII line in ( Oph fFig. ll2> the flux in 
the fragmented wind model is about 30% higher than that in the 
smooth wind model. Unfortunately, no models are available yet that 
would be able to predict the line fluxes based on detailed hydrody- 
namic simulations. A comparison of fluxes between modeled and 
observed lines would be a very sensitive tool. 



8.2 On tlie influence of the velocity field on the line profile 

The influence of the velocity field on model li ne profiles is il- 
lustrat ed in Fig. 1151 Hydrodynamic simulations of lFeldmeier et alj 
il997|) indicate that the velocity of the hot plasma is nearly the same 
as the cool-wind velocity. In this paper, we approximate the ve- 
locity field as monotonic ( see iFeldmeier & Nikutt j i2003) for the 
treatment of non-monotonic velocity). From the width of the line 
profiles we infer that the X-ray emission in all four stars is origi- 
nating in a region extending from about ~ 0.3i7oo to « O.Svoo- The 
corresponding radii of emission can be obtained if the velocity-law 
parameter /? is known. Vice versa, if radii of emission are given, /3 
can be inferred from the line-profiles. Usually, /3 is estimated from 
UV lin e fits, and typically lies between 0. 5 and 1.5. IFullerton et alj 
J2006h obtained for (Pup l3 = 0.5, while IPuls et al.N200el) obtain 
/3 = 0.95 and notice that different values may apply for inner and 
outer parts of the wi nd. /3 ~ 1 was inferre d from the line profile 
variability in CPup iEversberg et alj[l998t) . It seems that optical 
and UV line profiles are less sensitive to the exact value of /3 than 
the X-ray lines studied here. As shown in Fig. 1151 the line profile 
changes significantly in dependence on /3. The synthetic lines for 
our four program stars shown in the figures were actually slightly 
optimized by adjusting /? within the uncertainties of previous esti- 
mates (see figure captions). 



8.3 On the applicability of smooth wind radiative transfer in 
clumped winds 

The optical depth in a smooth wi nd is given by Eq . jsl, w hich was 
employed bv lKramer et aljj2003h and lCohen et aljj200d) to fit the 

line profiles observed in the X-ray spectra of ( Pup and C Ori. They 
have used tq = Af/47ruoo as one of their four free fitting param- 
eters. Assuming — 125 g cm-^ for all analyzed lines in Ori, it 
was shown that the best-fitting values of tq imply significantly re- 
duced mass-loss rates, e.g. due to the wind clumping. In a clumped 
wind, however, the optical depth is given by Eq. Jl 1> . which is more 
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Figure 15. Observed (histogram, co-added MEGitl) and modeled NeX 
line in Pup. From left to right line profiles are calculated assuming 
/3 = 0.7, 0.9, 1.5. 



general than Eq. Js}- The latter applies only for a limiting case of 
numerous and optically thin clumps. 

Little is known so far about mass, density and geometry of in- 
dividual clumps. Clump optical depths are likely to be different at 
different radii in the wind. This is reflected in our stochastic wind 
model, where the mass of a clump is conserved, but its radial optical 
depth scales as r~ (Eq.|8}. Furthermore, since clump optical depth 
is proportional to k^, a clump can be optically thick at one wave- 
length and optically thin at another. As can be seen from Tabled 
the mass absorption coefficient differs by one order of magnitude 
between the wavelenghts of the Si XIV and the O VIII line. Figure|8| 
shows model lines of Si XIV and O VIII in ( Pup. For the Si XIV line 
the difference between clumped and smooth model is small, indi- 
cating that the clumps are optically thin at A 6.18 (Si XIV). On the 
other hand, the model lines of O VIII for smooth and clumped wind 
are significantly different, and, therefore, at least some fraction of 
the clumps are optically thick at A 18.97 (O VIII). Consequently, 
Eq. <1 1> can be applied for the Si XIV line, but cannot be used for 
the O VIII fine. 

Overall, when stellar mass-loss rates are revised down because 
of the wind clumping, the clumped wind radiative transfer should 
be applied for consistency. 



8.4 On the line profiles and clump geometry 

As discussed in detail in lOskinova et alj i2004) . the optical depth 
in a clumped wind is determined by two factors. The first factor 
is the probability that a photon which encounters a clump is ab- 
sorbed by clump material. This probability is obviously equal to 
one for optically thick clumps, Td S> 1, and is reduced otherwise 
to 1 — exp(— Tci) (see Eq. I1 H . 

The second factor is the probability that a photon propagating 
along the ray towards the observer does encounter a clump. When 
clumps are distributed randomly in the wind, this probability de- 
pends on the length of the path. The longer a photon travels through 
the wind, the higher is the probability that it will be intercepted by a 
clump. Therefore, photons emitted in the part of the wind which is 
receding from the observer have higher probability to be absorbed. 
Thus there are less photons in the red wing of a line reaching an 
observer than in the blue part. This is the basic physical reason for 
the blueshift of the X-ray emission line profiles. 
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Figure 16. Two model lines, both calculated with no = 4.1 X 10~ s~ 
and TO = 10. Flux is normalized to the unabsorbed line flux. The only 
difference between the models is the assumption of isotropic opacity (i.e. 
clumps as "balls") and anisotropic opacity (i.e. clumps as "pancakes"). The 
solid vertical lines mark the positions of line maxima. Isotropic opacity 
leads to a more blueshifted and significantly skewed line profile. 



Moreover, the probability that a photon does encounter a 
clump depends also on the clump cross-sections. Approximating 
clumps as spheres, the clump cross-section is the same from all di- 
rections, which is conceptually similar to the atomic cross-sections. 
The emergent line profiles are skewed in both cases. 

To explain the observed near-symmet ry of line profiles in 
the framework of the smooth wind model, Ixramer et alj i2003ft 
and lCohen et alj|2006l) suggested the reduction of mass-loss rates, 
which implies lower densities. When the particle number density 
is reduced, the photon mean free path is increased. In a clumped 
wind with spherical clumps, each clump is analogous to a parti- 
cle. In order for the line profile to be symmetric, the "particle" (i.e. 
clump) number density must be reduced, which consequently leads 
to a larger photon mean free path (i.e. porosity length). 

This deep analogy between smooth winds on one side and 
clumped winds with spherical clumps on the other side ex- 
plains the princiml similarity between the conclusions reached in 
[Cphen et al. ( 2009) and inlowocki & Cohen 1 20061) , namely the re- 
quirement of large mean free paths in order to produce symmetric 
emission line profiles. 

A qualitatively different behavior of the radiative trans- 
fer is achieved, however, for anisotro pic opacity (i.e . non- 
sph erical, oriented clump s), as studied in lFeldmeier et alj ^2003^ 
and lOskinova et al1i2004) . This effect of clump geometry on emer- 
gent line profiles is demonstrated in Fig. ll6l Two model lines are 
calculated with the the same fragmentation frequency and the same 
value of To. The only difference is that in one case isotropic opac- 
ity, i.e. spherical "ball-like" clumps are assumed, while the other 
model is calculated for the anisotropic opacity in the limiting case 
of clumps which are flat like "pancakes". Figure^] clearly shows 
that in the latter case the line is non-skewed and less blueshifted 
than in the former case. 

The assumption of flat clumps was used in the present paper 
to model the observed line profiles. We have demonstrated that this 
model yields good agreement with the observed data. 
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9 CONCLUSIONS 

The effect of clumping on the wind attenuation is threefold. Firstly, 
it reduces the wind absorption compared to the smooth-wind 
model, allowing for more flux to escape. Secondly, it leads to a 
weaker dependence of the effective opacity on wavelength, making 
it entirely grey in the case of optically thick clumps. This prediction 
of similar opacity for different lines agrees well with the similarity 
of line profiles observed at different wavelengths. Thirdly, it af- 
fects the shape of line profiles, but only in the case when effective 
anisotropic opacity is assumed. The radially compressed, slab-like 
shell fragments lead to the emergence of the symmetric, but blue- 
shifted lines profiles which are indeed observed. 

Inhomogeneity of the stellar wind can reduce the mass-loss 
rates estimated empirically from the analysis of recombination and 
Hq lines. These reduced mass-loss rates should be used in conjunc- 
tion with the stochastic wind radiative transfer formalism in order 
to consistently model X-ray emission line profiles. 

The radii in the wind where continuum optical depth for the 
X-rays is unity is below 10 i?, in Chandra MEG energy band 
for Pup, C, Ori, ^ Per, and C, Oph. Stellar wind in C, Oph is almost 
transparent for the X-rays, consequently the emission line profiles 
show no effect of wind absorption. We estimate that 0.01% of 
the mechanical energy of the wind of C, Oph is converted into X- 
ray emission. The X-ray emission lines observed in the spectra of 
C, Pup, C, Ori, and ^ Per are affected by wind absorption. This result, 
as well as shape of line profiles, is in agreement with the analysis 
of line ratios in the He-like ions. The above analysis constrains the 
origin of X-ray emission to the acceleration zone in the stellar wind. 

Summarizing, we conclude that by using independently de- 
termined parameters as input for the model of an inhomogeneous 
stellar wind we are capable to reproduce the observed X-ray emis- 
sion lines in O stars. The wind clumping explains the shape of the 
line profiles, i.e. their width, symmetry and blueshift, as well as 
their similarity across the spectrum. 
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